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It is argued that a multicomponent structure of e + e~ and pp high 
energy collisions determines the shapes of multiplicity distributions 
and oscillations of H q moments. QCD and convoluted NBD are used 
for comparison with experimental data. Predictions at LHC and 
higher energies are shown. Difference between e + e~ and pp is dis- 
cussed. 

Multiplicity distributions are the integral characteristics of multiparticle 
production processes. They can be described either in terms of probabilities 
P(n, E) to create n particles at energy E or by the moments of these distri- 
butions. It has been found that their shapes possess some common features 
in all reactions studied. At comparatively low energies, these distributions 
are relatively narrow and have sub-Poissonian shapes. With energy increase, 
they widen and fit the Poisson distribution. At even higher energies, the 
shapes become super-Poissonian, i.e. their widths are larger than for Poisson 
distribution. The width increases with energy and, moreover, some shoulder- 
like substructures appear. 

Their origin is ascribed to multicomponent contents of the process. In 
QCD description of e + e~-processes these could be subjets formed inside 
quark and gluon jets (for the reviews see, e.g., PP, j2]). In phenomenological 
approaches, the multiplicity distribution in a single subjet is sometimes ap- 
proximated by the negative binomial distribution (NBD) first proposed for 
hadronic reactions in j^j. For hadron-initiated processes, these peculiarities 
are also explained by the multicomponent structure of the process. This is 
either multiladder exchange in the dual parton model [I], [5], varying number 
of clans P or multiparton interactions [Zj, |Hj, 

Such evolution of the multiplicity distributions can be quantitatively de- 
scribed by the energy behavior of their moments. To introduce them, let us 
write the generating function of the multiplicity distribution as 
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In what follows, we will use the so-called unnormalized factorial J- q and 
cumulant K, q moments defined according to the following formulas 



T q = ]T P(n)n(n - l)...(n - q + 1) = Z ' 



d q \uG{E,z) 
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They define correspondingly the total and genuine correlations among the 
particles produced (for more details see [TU], Since both T q and KL q 

strongly increase with their rank and energy, the ratio 

first introduced in [TT], is especially useful due to partial cancellation of these 
dependences. The factorial moments T q s are always positive by definition 
(Eq. (J2J)). The cumulant moments JC q s can change sign. They are equal to 
for Poisson distribution. 

The generating functions for quark and gluon jets satisfy definite equa- 
tions in perturbative QCD (see [12], [2]). It has been analytically predicted 
in gluodynamics JJJ that at asymptotically high energies H q moments are 
positive and decrease as q~ 2 . At Z° energy they should have a negative 
minimum at q w 5. Some hints to possible oscillations of H q vs q at higher 
ranks at LEP energies were obtained in ^JJ. The approximate solution of 
the gluodynamics equation for the generating function ^3] agrees with this 
and predicts the oscillating behavior at higher ranks. These oscillations were 
confirmed by experimental data for e + e~ and hadron-initiated processes first 
in (TJ], later in J3] and most recently in [T^]. The same conclusions were 
obtained from exact solution of equations for quark and gluon jets in the 
framework of fixed coupling QCD [17 . A recent exact numerical solution of 
the gluodynamics equation in a wide energy interval ^H] coincides with the 
qualitative features of multiplicity distributions described above. 

In parallel, the NBD-fits of multiplicity distributions were attempted [H], 
[19J. The single NBD is parameterized as 

T(n + k) (m\ n ( m\-' n - k , , 

p "* E '= fi^a (t) ( 1+ t) (5) 
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where T denotes the gamma-function. This distribution has two adjustable 
parameters m(E) and k(E) which depend on energy. However the simple 
fit by the formula (0) is valid till the shoulders appear. In that case, this 
formula is often replaced [S| by the hybrid NBD which simply sums up two or 
more expressions like (jSj). Each of them has its own parameters rrij, kj. These 
distributions are weighted with the energy dependent probability factors Wj 
which sum up to 1. Correspondingly, the number of adjustable parameters 
drastically increases. The analytical expressions for first five moments of 
two-NBD fits have been derived |2Uj . 

It was proposed recently [21] that hadron interactions can be represented 
by a set of Independent Pair Parton Interactions (IPPI-model) and each 
of binary parton collisions is described by NBD with the same parameters. 
Then their convolution leads to a common distribution 



Jmax 



P(n;m,k) = E WjPNBD(n;jm,jk). (6) 

This is the main equation of IPPI-model. One gets a sum of negative binomial 
distributions with shifted maxima and larger widths for a larger number of 
collisions. No new adjustable parameters appear in the distribution for j 
pairs of colliding partons. The probabilities Wj are determined by collision 
dynamics and, in principle, can be evaluated if some model is adopted (e.g., 
see j22], El)- O ne can expect that at very high energies Wj is a product of j 
probabilities W\ for one pair. Then from the normalization condition 

Jmax Jmax 

E *>j = E < = i (7) 

3=1 3=1 

one can find W\ if jmax, which is determined by the maximum number of 
parton interactions at a given energy, is known. In asymptotics, j max — ■» oo 
and W\ = 0.5. 

The factorial moments of the distribution © are 

T(jk + q) (my (my 

with 
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T(jk + q) 
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The m-dependence of cumulant moments is similar. For H q moments one 
gets 

Thus, according to Eq. (fTUj) H q are functions of the parameter k only and do 
not depend on m in IPPI-model. This remarkable property of H q moments 
provides an opportunity to fit them by one parameter. 

Once the parameter k is found from fits of H q , it is possible to get another 
parameter m rewriting Eq. (jHJ) as follows 
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m = k[^-\ . (11) 

This formula is a sensitive test for the whole approach because it states 
that the definite ratio of g-dependent functions to the power 1/q becomes 
g-independent if the model is correct. Moreover, this statement should be 
valid only for those values of k which are determined from H q fits. Therefore, 
it can be considered as a criterion of a proper choice of k and of the model 
validity, in general. 

We have compared IPPI-model conclusions with experimental multiplic- 
ity distributions of E735 [2H] collaboration for pp collisions at energies 300, 
546, 1000 and 1800 GeV extrapolated [21], |23 to the full phase space. An 
analysis of experimental data done in [24] has shown that 2 parton pairs are 
already active at energies above 120 GeV. We assume that 3 parton pairs are 
active at 300 and 546 GeV, 4 at 1000 and 1800 GeV, 5 at 14 TeV and 6 at 
100 TeV with NBD for a single collision. We divide by 2 the charged particle 
multiplicity and consider multiplicity of particles with the same charge (say, 
n+). Then, factorial and H q moments are obtained from experimental data 
on P(n+) according to Eqs. ©, ©• 

Experimental H q moments were fitted by Eq. ()10|) to get the parameters 
k(E) of the IPPI-model. We show in Fig. 1 how perfect are these fits at 
1.8 TeV for k equal to 3.7 and 4.4. With these values of the parameter k, 
we have checked whether m is constant as a function of q. Experimental 
factorial moments and IPPI- values for <p q were inserted in Eq. (II 1J) . The 
m(q) dependence is shown in Fig. 2 for the same values of k and for much 
larger value 7.5. The constancy of m « 12.94 is fulfilled with an accuracy 
better than 1.5% for k = 4.4 up to q = 16. 
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The same-charge multiplicity distribution at 1.8 TeV has been fitted with 
parameters m = 12.94 and k—4.4 as shown in Fig. 3. The similar procedure 
has been applied to data at energies 300, 546, 1000 GeV. We have found 
that m logarithmically increases with energy. The parameter k is equal to 
4.4 at 546 GeV and about 7-8 at 300 and 1000 GeV. Extrapolating m and 
choosing k=4A and 8, we calculated multiplicity distributions at 14 TeV 
and 100 TeV shown in Fig. 4. Also shown, is its shape at 14 TeV if ladder 
model probabilities [22], |H] are inserted in Oscillations of H q still persist 
|21j . Their amplitudes (~ 10 -2 ) are about two orders of magnitude larger 
than those for e+e~ at Z° (« 10^ 4 ). No suitable IPPI-fit by Eq. ((TDJ) has 
been found for e + e~ data. The asymptotical behavior of H q moments in 
IPPI- model is incompatible with 1/q 2 dependence predicted by QCD. The 
anomalous fractal dimensions of e + e~ and pp also differ ^D]. Thus, there is 
no direct similarity of e + e~ and pp (see [21] for more details). 

To conclude, the model of independent pair parton interactions is quite 
successful in describing experimental multiplicity distributions. Oscillations 
of H q moments can be ascribed to the mutlicomponent structure of both 
lepton- and hadron-initiated processes. Predictions for LHC energies are 
presented. No complete similarity of multiplicity distributions in e + e~ and 
pp has been observed. 

This work is supported in part by the RFBR grants 02-02-16779, 03-02- 
16134, NSH-1936.2003.2. 

Figure captions. 

Fig. 1 A comparison of H q moments derived from experimental data at 
1.8 TeV (squares) with their values calculated with parameter k=4A (dash- 
dotted line) and 3.7 (solid line). 

Fig. 2 The g-dependence of m for k=4A (squares), 3.7 (circles) and 7.5 
(triangles). 

Fig. 3 The same-charge multiplicity distribution at 1.8 TeV, its fit at 
m = 12.94, k = 4.4 (solid line). The dash-dotted line demonstrates what 
would happen if NBD is replaced by Poisson distribution. 

Fig. 4 The same-charge multiplicity distributions at 14 TeV and 100 TeV 
predicted by IPPI-model with logarithmic extrapolation of m (solid line - 14 
TeV, fc=4.4; dash-dotted - 14 TeV, k=8; dashed - 100 TeV, k=4A; for the 
ladder model: dotted - 14 TeV, k=4A). 
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